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Abstract. For the spherical model with nearest-neighbour interactions, the 
microcanonical entropy s{e,m) is computed analytically in the thermodynamic 
limit for all accessible values of the energy e and the magnetization m per spin. 
The entropy function is found to be concave (albeit not strictly concave) , implying 
that the microcanonical and the canonical ensembles are equivalent, despite 
the long-range nature of the spherical constraint the spins have to obey. Two 
transition lines are identified in the (e, m)-plane, separating a paramagnetic phase 
from a ferromagnetic and an antiferromagnetic one. The resulting microcanonical 
phase diagram is compared to the more familiar canonical one. 
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1. Introduction 

Ensemble averages are at the very core of equilibrium statistical mechanics as 
introduced by Boltzmann and Gibbs more than a century ago. The crucial idea is 
that long-time averages, which are the quantities of interest in equilibrium statistical 
physics, coincide with ensemble averages, with the advantage that the latter are easier 
to compute. An ensemble average is obtained by assigning statistical weights to all 
the microstates of a system, and calculating expectation values with respect to these 
weights. The choice of the statistical ensemble, i. e. the choice of statistical weights, 
depends on the physical situation that one is interested in: The microcanonical 
ensemble for example is appropriate for the description of an isolated system at 
fixed energy, whereas the canonical ensemble describes a system in equilibrium 
with an infinitely large heat bath of temperature T. For a suitable class of short- 
range interactions, both ensembles are known to give equivalent results in the 
thermodynamic limit (see [T] for details). For actual calculations, equivalence of 
ensembles allows one to use the ensemble which appears to be the most convenient, 
independently of whether it reflects the physical situation of interest or not. In this 
spirit, numerical calculations of microcanonical quantities have been used repeatedly 
to identify first order phase transitions, as the finite-size precursors of these transitions 
appear to be more pronounced in the microcanonical ensemble [H |21 HI [S] . Analytic 
calculations in the microcanonical ensemble, however, turn out to be more demanding 
than in the canonical ensemble, as Gibbs already remarked in his seminal treatise 



Microcanonical entropy of the spherical model with nearest-neighbour interactions 2 

"It is sufficient here to remark that analytically the canonical distribution is 

much more manageable than the microcanonical. " 
As a consequence, rather few such microcanonical calculations can be found in the 
literature. 

In the present paper, a microcanonical analysis is reported for the spherical 
model with nearest-neighbour interactions on a d-dimensional hypercubic lattice. For 
d ^ 3, this model is a toy model of a ferromagnet, showing a temperature-driven 
continuous phase transition from a ferromagnetic to a paramagnetic phase. Exposed to 
a magnetic field and at sufRciently low temperature, the model also undergoes a field- 
driven discontinuous phase transition. Although the nearest-neighbour interactions 
are clearly of short-range nature, the model is subject to a spherical constraint 
(from which it derived its name) which is effectively long-range. In contrast to the 
short-range interacting case, it cannot be taken for granted that microcanonical and 
canonical calculations yield equivalent results for long-range interacting systems. As 
a consequence, it was not clear to the author from the outset whether or not to 
expect equivalence of microcanonical and canonical calculations for the spherical model 
[71 [8] . The possibility to observe such nonequivalence renders the spherical model an 
interesting case for a study in the microcanonical ensemble. 

Behringer j9j has published a calculation of this model's microcanonical entropy 
s(e,m), considered as a function of the energy e and the magnetization m per spin. 
However, his analysis is valid only in a certain region of the (e, m)-plane, leaving 
out the parameter values in the phase-coexistence region where nonequivalence of 
ensembles might possibly occur. In this region, a saddle point analysis of the 
asymptotic integral which yields the microcanonical entropy is no longer possible, 
and it is one of the main results of the present article to extend the calculation of the 
microcanonical entropy into the coexistence region. 

Apart from completing the microcanonical calculation of s(e, m) for the spherical 
model, the results reported in the present article also serve a more general purpose: 
from the asymptotic analysis of a Laplace-type integral, one observes that the 
occurrence of a phase transition is related to the breakdown of the above mentioned 
saddle point analysis. A careful study of this asymptotic integral nicely illustrates 
the mathematics at work which leads to a nonanalytic behaviour of the entropy, and 
hence to a phase transition. 

Finally, from the calculation of the microcanonical entropy s(e, m) for all 
accessible values of e and m, we obtain the microcanonical phase diagram of the 
spherical model, i. e. the lines in the (e, m)-plane for which s is not a smooth function 
and phase transitions occur. This diagram comprises, as expected, a transition line 
corresponding to the ferromagnetic to paramagnetic transition of the spherical model, 
but it also contains a second transition line, signalling a transition from a paramagnetic 
to an antiferromagnetic phase. Omitting many of the details, this phase diagram, and 
a similar one for the Ising model, have been published previously as a Letter |10) . 
suggesting that the observed behaviour is typical for short-range ferromagnets in the 
microcanonical ensemble. 

The rest of the article is organized as follows. In section [2l the spherical model 
is introduced, and a calculation of the microcanonical entropy s{e, m) of this model is 
reported in section[3l The resulting caloric curves and response functions are discussed 
in section|4l From the results for the microcanonical entropy, the canonical free energy 
and the canonical phase diagram are recovered by means of a Legendre transform in 
section [5j Finally, a summary and a discussion of the results is given in section [G] 
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2. Spherical model with neairest-neighbour interactions 

Consider a hypercubic subset 

{1,...,LY d^'^ (1) 

of a d-dimensional hypercubic lattice. On each of the N = L"^ lattice sites a real degree 
of freedom (T^ e R is placed. The spherical model is characterized by the Hamiltonian 

i?Ar : R^ ^ R, -J^aiaj, (2) 

where J > is a coupling constant and the angular brackets denote a summation over 
pairs of nearest neighbours on the lattice C. This model was introduced by Berlin 
and Kac [TT] as an exactly solvable caricature of the Ising model of a ferromagnet. 
In contrast to the Ising model case, the "spin variables" di are real numbers: their 
modulus is not fixed to unity as for the Ising model where Ci € {— 1,+1}. Instead, 
the spherical constraint 

N 

Y.^^ = N (3) 

1=1 

allows for fluctuations of the modulus of the spin variables. In the canonical ensemble, 
this model has been solved in the thermodynamic limit for arbitrary d, and a 
continuous phase transition from a ferromagnetic phase at low temperatures to a 
paramagnetic phase at high temperatures (or energies) occurs for all d 5^ 3 '12]. 



3. Microcanonical entropy 

For the spherical model with nearest-neighbour coupling, an analytic calculation of 
the microcanonical entropy as a function of energy e and magnetization m has been 
reported by Behringer in [9] , but the analysis does not apply to certain regions in the 
(e, m)-plane. Starting point for such a microcanonical calculation is the density of 
states 

OA,(e,m)= / da5[Ne--HN{(j)]5[Nm-MN{(j)]5(N-y"ai\,{'i) 
where the function 

N 

Mat : R^ R, ^ XI 

i=l 

yields the magnetization of a microstate. It was shown in [S^ that, asymptotically for 
a large number N of lattice sites, the density of states can be written as 

/•a+ioo 1 /"C+icxD i /•6+ioo i 

^ , , / az I au I aw .^^ , , 
VtNie.m)^ — — — exp[7V(?!)^,„,„(£,m)],(6) 

J a—ioo ^'^ J c—ioo ^'^ J b—ioo '^'^ 



where 

(e, m) = ze + wm + u + — — — — — - ^ / ,^ f , \n i u - zjS^ cos ipj ] .(7) 
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From this expression, the microcanonical entropy in the thermodynamic hmitjj 

s{e,m) = lira — hi iljv(e, m), (8) 

can be computed by the method of steepest descent (see section 6.6 of reference 
[13]). Bchringer proposes to do that by studying the saddle points of (t)z,w,u in the 
complex {z,w,u)-space (see equations (28)-(30) of reference [9]), and for his purposes 
such an analysis turns out to be sufficient. To obtain the entropy s in the entire 
range of accessible values in the (e, r7i)-plane, however, a saddle point analysis is not 
satisfactory. The reason for that, as we will see in the following, is that branch cuts 
in the complex (z,M)-space due to the logarithm in Q come into play. 

We start by observing that the expression for the density of states can be 
simplified by performing the w- and the M-integration in ^ , yielding 



J7Ar(e, m 



X exp - — 7 / dV In , ;~ ^ (9) 

(the derivation of this result is given in [Appendix Ap . Due to the logarithm in 
the integrand of the z-integration has two branch cuts on the real line, one extending 
from — oo to —l/{dJ), the other from +l/{dJ) to +oo. Apart from these branch cuts, 
the integrand is holomorphic, and the contour of integration can be deformed freely, 
as long as it does not cross the cuts. For an asymptotic evaluation of the integral in 
the large- limit by means of the method of steepest descent, the path of integration 
in the complex plane is deformed such that its imaginary part becomes zero. In this 
limit, the value of the integral is given by the integrand of the z-integration in ^ 
evaluated at the maximum along that path. This maximum is located in the interval 

I =[-l/idJ),+l/idJ)] (10) 

of the real z-axis, and — depending on the values of s and m — it may be either a saddle 
point of the exponent in or one of the boundary points ±l/(c?J). 

For lattice dimensions d — I and d = 2, the maximum of the integrand along the 
path of vanishing imaginary part is always located at a saddle point in the interior 
of X. As we will see from the following discussion, such a behaviour corresponds to 
the absence of a phase transition. For d ^ 3, however, the situation is different: In 
figure [U the behaviour of the paths of vanishing imaginary part upon variation of 
e and m is illustrated for the spherical model on a lattice of dimension d — 3 with 
coupling constant J = 1, but the behaviour is similar for larger values of d. We can 
identify three different regions in the (e, TO)-plane which are sketched in figured 

Region I: For relatively small energies e, the maximum of the integrand of the 
z-integration is located at +l/{dJ), and the microcanonical entropy in the 
thermodynamic limit is given by 

sHe, r.) 4 m [2e(l + ^) ] - ^ /^^^^ In (l - i g cos . (11) 

Physically, region I corresponds to the coexistence region in which domains of 
positive magnetization and domains of negative magnetization coexist. The 

I Note that here and in the following Boltzmann's constant is set to unity. 
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Figure 1. Paths of vanishing imaginary part of the exponent of ^ in the complex 
2-plane. The paths shown here are for the spherical model on a three-dimensional 
lattice with coupling constant J = 1 and magnetization m = 0, and for various 
values of the energy e as indicated in the plot. One can observe in this figure 
that, for intermediate values of s, the paths cross the real axis in the interior of 
the interval [—1/3,-1-1/3]. For larger or smaller energies, the paths "get stuck" at 
the branch cuts at —1/3 and -1-1/3, respectively. 




Figure 2. Microcanonical phase diagram of the spherical model on a three- 
dimensional cubic lattice. The microcanonical entropy s is defined on the grey 
shaded region in the (e, m)-plane. Within each of the regions I, II, and III, the 
entropy is analytic, but not so on the boundaries separating the regions. 



independence of on the magnetization m in equation (jlip is a hallmark of this 
behaviour. This m-independence can be seen as a consequence of a remarkable 
interplay of numerator and denominator of the argument of the logarithm in ([9]) : 
The numerator determines the locations of the branch cuts (— oo, — l/((iJ)] and 
[+l/{dJ),+oo) in the complex z-plane, and precisely at the endpoint +l/{dJ) 
of one of these cuts the denominator, and therefore the entire exponent in ([5]), 
becomes independent of the magnetization m. 

Region II: For intermediate energies, the integrand in ^ has a saddle point zq in 
the interior of the interval T, determined by the equation 

d°(p 







+ m^dJ + [l - + zo(e + rn^dJ)] 



1 - zoJJ2-i^i cos ipj 



(12) 
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The microcanonical entropy in the thermodynamic limit is then given by 
""(e,™) = ^ In 2c(^l - TO^ + zo{e + TO^dJ)^j 



s 



i=i 

Region III: For larger energies e, the maximum of the integrand of the z-inte- 
gration in Q is located at — l/((iJ), and the microcanonical entropy in the 
thermodynamic limit is given by equation (jl3p with zo = — 1/ {dJ) . 

The boundaries between the different regions in the (e, m)-plane are given by the 
values of (e, m) for which the saddle point of the integrand in ([9]) approaches one of 
the boundary points zq = ±l/((iJ). These points are obtained by solving equation 
(|12p with zq = ±l/((iJ), yielding for the two boundary curves the parabolas 

m2 + ad[m2±(l-m2)] 
e±(m) = -dJ ^ (14) 

t + Od 

with 

J[<d,T:Y ^ « - Z^j=l COS </3j 

Since the microcanonical entropy s{e,m) in the thermodynamic limit is defined 
piecewise in the regions I, II, III, one expects this function to be nonanalytic (i. e. not 
infinitely-many times differentiable) at the borders of the regions. As a consequence, 
the plot in figure [2] can be interpreted as a microcanonical phase diagram, showing 
at which values of the control parameters e and m thermodynamic singularities are 
encountered. From an overall plot of the graph of s{e,m) as shown in figure [3] this 
nonanalytic behaviour is not immediately visible and we will consider derivatives of s 
to illustrate these features in the next section. 

In summary, we have observed that the asymptotic evaluation of the integral 
governing the density of states yields different behaviours for the regions 1, II, and 111, 
corresponding to the different phases of the spherical model. Only for region II is a 
saddle point analysis of the integral in ^ applicable. Upon reaching the boundaries 
of region II, the saddle point on the real z-axis gets stuck at the branch cut of the 
exponent in equation ([9|), leading to nonanalytic behaviour of the microcanonical 
entropy s(e, m). Remarkably, at the boundary e+{m) between regions I and II, the 
saddle point gets stuck at a point precisely such that the m dependence of the exponent 
in Q cancels, leading to an affinc (in the m-direction) region of the entropy. 



4. Caloric curves and specific heat at fixed magnetization 

Caloric curves or response functions like the specific heat are typical quantities which 
are measured in experiments, and kinks, discontinuities, or singularities in these 
quantities are hallmarks of a phase transitions. From the microcanonical entropy, 
these quantities are obtained as first or second derivatives (or functions thereof). To 
illustrate the occurrence of nonanalyticities in the microcanonical entropy s{e,m), 
we consider the spherical model at fixed magnetization toq- The corresponding 
microcanonical entropy function is 

Smois:) = s(e,mo), (16) 
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Figure 3. Graph of the microcanonical entropy s{e, m) of the spherical model 
with coupling strength J = 1 on a three-dimensional cubic lattice in the 
thermodynamic limit. The constant (in the m-direction) part on the left hand 
side of the graph corresponds to region I where phase coexistence occurs. 




Figure 4. Microcanonical caloric curves of the spherical model on a three-dimen- 
sional cubic lattice with fixed magnetization mo = 0.92 in the thermodynamic 
limit. Left: Inverse temperature /3 as a function of the energy e as defined in 
II17I I. Right: A more familiar representation of the caloric curve, showing e as a 
function of the temperature T = 1//3, restricted to positive values of T. 



and the microcanonical caloric curve is given as 

Pis) = ^^^, (17) 

where (3 — 1/T is the inverse temperature. For the magnetization toq = 0-92 the 
caloric curve is plotted in figure ID The nonanalyticity of the microcanonical entropy 
s is clearly visible only in the second derivative 

= (18) 
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Figure 5. Left: Second derivative s"ip{e) of the entropy at fixed magnetization 
mo = 0.92 for tlie spherical model on a three-dimensional cubic lattice in the 
thermodynamic limit. Right: Microcanonical specific heat c(e) as defined in llf 9|l 
for the same model and parameter values. 



or in the microcanonical specific heat 

both plotted in figure [5l Qualitatively, graphs for other values of the fixed 
magnetization mo look similar. 

The observed nonanalyticities of s(e,m) can also be interpreted as phase 
transitions in the spherical model with fixed magnetization. For any value of the 
fixed magnetization mo, two continuous phase transitions occur at the values e±(mo) 
as given in equation (fH)) . The order of the phase transition, i.e. the fact that it is 
a continuous one, could have been deduced also without computing the exact result 
for s{e,m) from geometric arguments similar to those applied to the Ising model in 
reference j 14) . 

5. Legendre transform and canonical free energy 

The microcanonical entropy that we have computed in section |2] forms the starting 
point of an analysis of the spherical model in the microcanonical ensemble. An 
analogous role is played by the canonical Gibbs free energy 

N 

9N{P,h)^-—l da<5(7V- Va2)e-^«"(-) (20) 



for calculations in the canonical ensemble. In the thermodynamic limit, the 
corresponding infinite-system quantity 

g(/3,/i)= lim gNW,h) (21) 

A/— >oo 

is related to the microcanonical entropy s(e, m) by means of a Legendre transform, 

- (3g{(3, h) = s{ei(3, h),m{(3, h)) - f3e{l3, h) + phrh{P, h), (22) 
where e(/3, h) and m(/3, h) are the solutions of the equations 

^s(e,m) -/3 = 0, ^s(e, m) -f /3/i = 0. (23) 
oe dm 

In this way we can derive the canonical result from the microcanonical one. While, in 

the thermodynamic limit, the canonical free energy is always given as the Legendre 
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Figure 6. Canonical phase diagram of the spherical model. The plotted lines 
in the /i)-plane mark the parameter values at which the canonical free energy 
h) is nonanalytic and phase transitions take place. The right line corresponds 
to the transition line e+ (m) in the microcanonical case at which a transition from 
paramagnetic to ferromagnetic behaviour takes place. The left line corresponds 
to the transition line £_(m) in the microcanonical case at which a transition from 
antiferromagnetic to paramagnetic behaviour occurs. 



transform of the microcanonical entropy, the inverse is not always true: Equivalence 
of ensembles holds only for concave entropy functions [El [16] , and only in this case 
can the entropy be recovered from the free energy by means of a Legendre transform. 

The equations in (|23p provide a way to translate the microcanonical control 
parameters (e, m) into the canonical control parameters {f3,h). We can use these 
relations to translate the microcanonical transition lines p4)) in the (e, m)-plane into 
the canonical transition lines in the /i)-plane. From the microcanonical transition 
line e+(m) in equation (|23p we obtain the canonical transition line /i+(/3) = for 
/3 G [(1 + ad)/2dJ, +oo), with as defined in equation p5|) . From the transition line 
e_(m) we obtain the canonical counterpart 

_ 8dJ{l + ad) 

^-^''^ ~ - {Adjy ^ ^ 

for ft, e R. Details of the calculation are given in [Appendix B[ and the resulting 
canonical phase diagram is plotted in figure [SI 

The two transition lines in the (/3, /i)-plane display the typical behaviour of a 
ferromagnetic transition (for positive /3) and of an antiferromagnetic transition (for 
negative and this observation deserves some remarks. For physical reasons, the 
canonical ensemble is often restricted to positive temperatures, corresponding to 
positive values of the inverse temperature /3. First, this appears to be appropriate 
for modelling a system in a thermal bath, but there is a second reason for this 
restriction: For many physical systems, in particular those with a standard kinetic 
energy term which is a quadratic form in the momenta, no canonical equilibrium 
states corresponding to negative temperatures exist. This is different for the spherical 
model, as well as for many other spin models for which the energy per degree of 
freedom is bounded above. 

Imagining the physical situation of a macroscopically large spherical model 
coupled to a thermal reservoir with positive inverse temperature only the right 
half-plane in figure[6]can be explored and the ferromagnetic-to-paramagnetic transition 
can be probed. The left half-plane, corresponding to negative temperatures, might 
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still be accessed by means of a trick: Since, in the canonical ensemble, only the 
product (3Hn of f3 times the Hamiltonian function H]\f enters the calculation of 
thermodynamic quantities, the thermodynamics of a system with Hamiltonian Hj^ 
at negative temperatures is identical to that of a system with Hamiltonian —Hj^ at 
positive temperatures. For the spherical model, this simply corresponds to an inversion 
of the sign of the coupling constant J in the Hamiltonian function ((2), giving rise 
to a spin system with antifcrromagnctic coupling J < 0. This reasoning explains 
the appearance of a transition line of antifcrromagnctic type in the left half-plane of 
figure O 

It also follows from this discussion that equivalence of the microcanonical and 
the canonical results for the spherical model requires the canonical ensemble to be 
considered for positive as well as negative temperatures: Only then it is possible 
to derive from the canonical free energy, by means of a Legendre transform, the 
microcanonical entropy s{e,m) on its entire domain, implying equivalence of the 
microcanonical and the canonical ensemble on the thermodynamic level [111 [12] • Note 
that this is a common situation for spin models, occurring for the Ising model, the 
Potts model, and many others. 

6. Discussion and conclusions 

We have presented a calculation of the microcanonical entropy s(e, m) of the spherical 
model with nearest-neighbour interaction, valid for all accessible values of the energy e 
and the magnetization m per spin. The entropy function is found to be concave (albeit 
not strictly concave) and, as a consequence, equivalence of the microcanonical and the 
canonical ensemble holds. Since, in addition to the nearest-neighbour interactions, the 
degrees of freedom of the spherical model are subject to an — effectively long-range — 
constraint, equivalence of ensembles could not have been taken for granted from the 
outset. We can interpret our finding of a concave entropy, concluding that the long- 
range constraint is not sufficient to change the overall characteristics of the model into 
that of a long-range interacting one. Only when the interactions between the spins are 
made properly long-range (like in the mean- field spherical model discussed in [T71ll8|) 
would we expect the entropy s(£, m) to develop a nonconcave part and noncquivalence 
of ensembles to occur. 

When computing the microcanonical entropy of the spherical model by means 
of the method of steepest descent in section |31 we have observed different regions, 
depending on the values of the parameters e and m. In region II, the main contribution 
to the relevant asymptotic integral comes from a proper saddle point on the real axis 
of the integration variable z. In regions 1 and 111, however, the saddle point gets 
stuck at a branch cut in the complex z-plane. It is this "getting-stuck" which leads to 
nonanalytic behaviour of s, giving rise to a phase transition in the spherical model. 

Note that a similar mechanism of generating nonanalytic behaviour from the 
interplay of a saddle point and a branch cut has been observed previously. In the very 
different context of nonequilibrium stationary states described by a Langevin equation, 
Farago [19] calculates the probability distribution function of the injection of energy 
into the system. His calculation also makes use of the method of steepest descent 
and, depending on some parameter values, a saddle point can approach a branch cut 
in the complex plane (Figure 3 of [H]). Once the saddle points gets stuck there, 
the probability distribution function develops a nonanalyticity, similar to what we 
observed for the microcanonical entropy of the spherical model. Farago also remarks 
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that this nonanalyticity can be interpreted as a phase transition; for details see the 
discussion in section 2.3 of reference [12]. 

The hnes in the (e, m)-plane at which the microcanonical entropy is nonanalytic, 
plotted in figure [H can be interpreted as a microcanonical phase diagram. Compared 
to its canonical counterpart shown in figure [H the microcanonical phase diagram looks 
remarkably different, and also its physical implications are unfamiliar when compared 
to the canonical situation: Varying, for example, the energy e while keeping fixed the 
magnetization at any value of m, one typically crosses two transition lines in the phase 
diagram, therefore observing two phase transitions, signaled by kinks in the specific 
heat. Similarly, four transition lines are crossed upon variation of the magnetization 
while keeping the energy e fixed at any value — d < e < £+(0), while two transition 
lines are crossed for energies e+(0) < £ < e_(0). Qualitatively similar results can also 
obtained for the Ising model |10j . and presumably for many other short-range spin 
models. For more details, especially on physical realizations of such models under 
microcanonical conditions, the reader is referred to reference |10) . 
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Appendix A. Derivation of equation {[9]) 

Starting from equation Q with exponent ([7]) , the w-integration is a Gaussian integral 
which is straightforward to perform: 



r2Ar(£,m) 



X exp I N 

a+ioo 



a+ioo 



dz 



c+ioo 



dw 
2^ 



w 



4(u - zdJ) 



d'^ip 

wm + ze + u— - I ,„ ■ , 



dz 



c+ioo 



du 



2^ / . 2^^'^''^ 

a~ioo ^" J c-ioo ^" 



N{fiu,z)-- 



6-f ioo 



dw 



a+ioo 



dz 



c+ioo 



1 



2vNTTJa-ioo 27r 7c-ioo 27r y^a{u,z) 
with 

1 



exp 



ioo 

N[f{u,z) 



27r 



exp 



4a(M, z) 







m ^ 






)] 


(A.l) 



a(u, z) 



4(u - zdJ) 



and 



1 /■ d'^(i5 / 

/(u, z) = u + ze- - \ ——g In u - z J V ( 

2 J[0,27r)<i K^T^r \ 



(A.2) 



(A.3) 



The Gaussian integral in (jA.ip has been solved assuming that 5R(m) > '3i{zdJ) along the 
contour of integration. To solve the u-integral in (jA.ip , we perform the substitution 
of variables 



xu, 



(A.4) 
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^ , , i f dx , — ( N f , / \^ M 

SZjv(e,m) ~ , / — V 1 — xdJexp < / - — -7 In 1 — xJ > coses,- > 



du 



-u 



{3-N)/2 [i + xe- m^{l - xdJ)] } . (A.5) 



2tt 

In a second step we substitute 

Nu[l + xe ~ m'^{l - xdJ)] y (A. 6) 

in the second integral in (|A.5|) . yielding 



iVvr J 27rY [l + a;e-m2(l-xdJ)] 



Af-5 



X exp 
where 



9(W) = / gs<»-"l'»e». (A.8) 

The y-integration is along some path extending from a — ioo to a + ioo for a > 0. 

For (N — 3)/2 e N, the integrand in (lA.Sp has a single pole at the origin, and by 
means of the residue theorem we obtain 

^^^^ " [{N-5)/2]l " r[(iV-3)/2]' ^^'^^ 

For (A^ — 4) /2 G N, the integrand in (jA.Sp has a branch cut along the negative real 
axis. The contour of integration is deformed such that — for some small S > — it goes 
from a — ioo to —00 — iS, follows slightly below the branch cut and back to —00 + iS 
slightly above the cut, and then to a + ic». Only the contours along the branch 
cut contribute to the integral, yielding integral representations of the F-function. In 
summary, we obtain again 

g(N) ^ —r\(5^ N)/2] = — ^——r (A.IO) 

' 2tt " ^ F[(Ar-3)/2] ^ ' 

Inserting these results into (|A.7[) and rearrangement of the terms leads to equation 



Appendix B. Derivation of the canonical phase diagram 

In this appendix we derive the canonical phase diagram as shown in figure [6] from 
the microcanonical one in figure [H For this purpose, we translate the microcanonical 
transition lines e± (to) from equation (|14p into the corresponding values of /3 and h by 
means of the equations in . 
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Transition line £+(m); This transition line separates the regions I and II as defined 
in section [31 On the boundary we can therefore use either s^ from equation pT|) or 
s^^ from equation (jl3p . and for simplicity we choose the former. Inserting s^ into the 
equations in we obtain 

(9 

— s^{e,m)=0^-Ph, (B.l) 
dm 

§-/ie.m) =^^^-m- (B.2) 

As a consequence, the canonical transition line is located at /i = 0. Evaluating equation 
(|B.2p at the microcanonical transition line (m) , we notice that e+ can take on values 
ranging from —dJ to e+(0) = —dJad/{l + a^), which implies that 

(3 e [f3{e+m, (3{-dJ)) = [i±^, +CX3) (B.3) 

with ad as defined in equation psp . 



Transition line £_ (m) ; This transition line separates the regions II and III. We choose 
to insert s^^^ — given by equation with zq = —l/{dJ) — into the equations in ([^5]) . 
yielding 

•») - - 2.,,(l-t'^)-2e - •»>• <«^*' 

= - 2A/(l-L^)-2e = «^-™)- (''■'5> 
Dividing (|B.4p by (jB.Sp we obtain 

/i = -AdJm. (B.6) 
Now we evaluate (jB.sp at e_(TO), 

and by inserting expression (|B.6p into (|B.7p we obtain the canonical transition line 
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